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Abstract 

A new renormalization group treatment is proposed for the critical exponents 
of an m-fold Lifshitz point. The anisotropic cases (m 7^ 8) are described by 
two independent fixed points associated to two independent momentum flow 
along the quadratic and quartic directions, respectively. The isotropic case is 
described separately. In that case, the fixed point is due to renormalization 
group transformations along the quartic directions. The new scaling laws are 
derived for both cases and generalize the ones previously reported. 
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Since the early formulation of the m-axial Lifshitz critical behavior associated to the 
onset of helical order in magnetic systems |1],|2[ , there have been various studies pointing out 
other applications of this kind of critical behavior in other real physical systems like high- 
T c superconductors 0, ferroelectric liquid crystals [|J, magnetic materials and alloys 
etc. These multicritical points appear at the confluence of a disordered phase, a uniformly 
ordered phase and a modulated ordered phase The modulated phase is characterized 

by a fixed equilibrium wavevector which goes continuously to zero as the system approaches 
the Lifshitz point. In case this wavevector has m-components, the critical system presents 
an m-fold Lifshitz critical behavior. If m — d, it displays an isotropic critical behavior. 
There are several types of isotropic behavior, but we will be concerned here only with 
the case m = d near to 8. Otherwise (m ^ d ^ 8), the system presents an anisotropic 
critical behavior. The particular interest in the isotropic case m = d near to 8 is due to its 
similarity to the anisotropic case, for both ones can be studied using the same field-theoretic 
description. 

The purpose of this letter is to present new renormalization group (RG) arguments and 
new scaling laws that point towards a natural explanation of the m-fold Lifshitz critical 
behavior for the anisotropic and isotropic cases. Our reasoning will be based on mag- 
netic systems. The m-fold Lifshitz point can be described by a generalization of the ANNNI 
model ||, which is a spin-| Ising model on a <i-dimensional lattice with nearest-neighbors in- 
teracting ferromagnetically as well next-nearest-neighbors with antiferromagnetic couplings 
along m directions. Its field theoretical representation is given in terms of a modified (ft 4 the- 
ory with higher derivative terms along the m directions, with the following bare Lagrangian 
density 

L =\\V 2 m<t>o\ 2 + \\ V(d-m) 00 I' + 5 ^| Vm 00 ? + ^0$ + ^0&t (1) 

At the Lifshitz point 5q = to = 0. The renormalized correlation functions are defined at 

-i 

the critical theory for t = 0, M — 0, where t = Z7 2 to, M = 0o are the renormalized 
reduced temperature and order parameter (magnetization), respectively . Below the Lifshitz 
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critical temperature Tl, one can expand the renormalized vertex parts for t ^ 0, M ^ 
around the ones for t = 0,M = 0asa power series in t and M, whereas above Tl the 
renormalized vertices for t 7^ are expanded around the renormalized parts calculated at 
t = in powers of t 0. We shall keep do = from now on, since this choice for the 
renormalized theories simplifies the dimensional analysis. The dependence on the external 
momenta along the quadratic (noncompeting) and quartic (competing) directions can be 
split in two independent contributions. The RG flow in the quartic momenta scale can be 



done separately, following a situation described by Wilson in the early seventies flO|| . The 
absence of quadratic terms in the components of the momenta along the competition axes 
allows us to perform a dimensional redefinition of the m-dimensional momenta subspace. 
This is the main advantage of choosing 5q = as the starting point for our RG analysis. 

Our new point is to consider the momenta along the competing directions entering with 
half the power of a momentum scale in dimensional analysis. For the anisotropic case, the 
momenta scales associated to the two orthogonal subspaces, namely the (d — m)-dimensional 
noncompeting one, characterized by the renormalization scale K\, and the competing m 
directions chacterized by k 2 , flow independently under RG transformations. This produces 
two independent fixed points, which lead to two new independent sets of scaling laws for 
each subspace. These new relations are valid to arbitrary order in an e^-expansion, where 
€l = 4 + y — d. In particular, the conventional Josephson hyperscaling relation for the 
anisotropic m-fold Lifshitz critical behavior [T|], 

2 - a L = (d - m)v L 2 + mv Li , (2) 

should be replaced by two independent new relations, namely the one associated to correla- 
tions perpendicular to the competition axes 

nrn 

2 - a L2 = (d - y)^2, (3) 
together with the hyperscaling law associated to correlations along the competing directions 

Tfl 

2-a L4 = 2{d- —)v M . (4) 
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Furthermore, the other new scaling relations appropriate for the critical exponents per- 
pendicular to the competition axes, 

7L2 = ^L2(2 - 771,2), (5a) 
x (d-f)+2- VL2 

h2 = RRT^ ' (5b) 

Pl2 = -v L 2{{d--)-2 + r ]L2 ), (5c) 

imply the Widom ^ L2 = Pl2^L2 — 1) and Rushbrook a i2 + 2Pl2 + 7L2 = 2 relations. The 
ones associated to critical correlations along the competition axes are analogous, namely 

7L4 = ^L4(4 - rj L4 ), (6a) 



2(d-f)+4-r ?L4 



(6b) 



2(d-f)-4 + 77 L4 ' 

1 Til 

/?L4 = -^4(2(d-y)-4 + 77 L4 ), (6C) 

which lead to 7^4 = Pla{^la — 1) and + 2/5^4 + 7^4 = 2. Note that at one-loop level, 
^L4 = \ v L2i and the above equations reduce to the scaling relations given in Ref. 0. In 
addition, for the isotropic behavior m = d near to 8 (the expansion parameter is ei = 8 — d) 
and the new corresponding scaling relations are: 

2 - a LA = dv Li , (7a) 

7L4 = ^L4(4 - rju), (7b) 
d + 4 - rj LA 



J LA 



d - 4 + 77x4 ' 



(7c) 



/?L4 = 2^4(^-4 + 77x4), (7d) 

which give rise to the same Widom and Rushbrook scaling relations mentioned above. 

Let us describe the method in some detail. We start with the anisotropic case. We 
choose two independent sets of normalization conditions for the critical theory, i.e., two 
different (subtraction) symmetry points. The first one is chosen with nonvanishing external 
momenta along the noncompeting (d— m)-dimensional subspace. The theory is renormalized 
at (quadratic) external momenta scale K\. The flow in this scale gives origin to the critical 



indices vl2 and r]L2- The second one is denned at nonvanishing external momenta along 
the m-dimensional competing (quartic) subspace. The theory is renormalized at (quartic) 
external momenta scale K2, originating the critical exponents and 77^4. Therefore, there 
are two sets of renormalized vertex parts, characterized by the scales K\ and k 2 , defined by: 

where g T are the renormalized coupling constants, Z^ and Z^^ are the field and tem- 
perature normalization constants and A T are the associated momentum cutoffs along the 
quadratic and quartic directions, respectively. (Except for which cannot be renormal- 

ized multiplicatively, all other vertex parts can be managed within this argument.) The label 
t = 1(2) refers to the nonvanishing quadratic (quartic) external momenta and zero quartic 
(quadratic) external momenta. As usual, ki (i = 1, N) are the momenta associated to the 
vertex functions T^' L ^ with N external legs and pi (i = 1, L) are the momenta associated 
to the L insertions of 2 operators. In terms of dimensionless parameters, g T = u t (kI t )~^ , 
and A r = uo t (k^. t )~ . As usual, uo T , Z^ and Z^^ can be represented as power series in 
u T . (Henceafter, we shall suppress the cutoff A T .) Invariance of the bare functions under the 
momenta scales used to fix the renormalized theories leads to the RG equation for T^ffi 

d 1 e 

where (3 T = («r§^), l<p(r){u T ) = Pr dln ^ T) and 7^ (t) (m t ) = -ft /'"^^ are calcu- 
lated at fixed bare coupling A T . In terms of the bare (dimensionless) coupling constant, 
f3 T = — re^( 8; g" 0T )~ 1 . The volume element in momenta space associated to loop integrals, 
d d ~ m qd m k, has engineering dimension A d_m A^ = A d ~^, where q([q\ = A) represents a 
(d — m) -dimensional vector perpendicular to the competing directions and k([k] — Aa) an 
m-dimensional vector along the competing axes, respectively. Under a flow in the exter- 
nal momenta k^ T ), we have the following simple scaling properties at the fixed points u* 
((7V,L)^(0,2)): 
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N(d— m ) tVyAf \ 

r (N ' L) (ok-, n OD; ML* k)~ AN+(d~f)~^j^-2L] ^ +Ll * {N>L) 

R(t) \P K i(r)} PPl(r)i k t) — P L R(t) \ K i(t)i Pi{t)i " T 5 K t)- 

(10) 

Now take L = 0. Above the Lifshitz critical temperature, the RG equation takes the 
form: 

1 

^o^ + ^a^r ~ 2 Nl<t>{T){Ur) + t^m^^S?) = °- 

At the fixed point (3 T {u* T ) = 0, and the solution can be written as: 

rg^fc^tx,^) = K7^F^\k i{T)j K T f^). (12) 

Define r\ T = 7^( r ), # r = — 7^a( r )> sucn that = VL2, V2 = ?7z4> an d so on. Under a flow in 
the external momenta, we have: 

i^VV)>~V)(p- 2T t)^). (13) 

By choosing p = ^(-4^) «t+2t ) an d replacing back in (13), the vertex function depends 
only on the combination k^ apart from a power of t. Therefore, we can identify the critical 
exponents v x = u L2 and u 2 = z/ i4 as 

v; 1 = 2t + 6 t . (14) 

Hence, at the fixed point all correlation functions (not including composite operators) 
scale at T > Tl, since they are functions of k^ only. For N = 2 we choose p = fcj( T ), the 
external momenta. As fcj( T ) — > 0, = x -1 and we can identify the susceptibility critical 
exponent: 

7 r = z/ r (2r - rj T ), (15) 

which are Eqs. ( [5a] ) and (pa]). Below T^, the renormalized equation of state at the fixed 
point is expanded in powers of M, which under a flow in the momenta turns out to be 
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We choose p to be a power of M: 

M ... 2 

and from the scaling form of the equation of state H(t, M) = M 5t f(— V), Eqs. (|5T5|), (|5c|) ; 
follow in a straightforward manner. 
The specific heat exponents can be obtained by analysing the RG equation for 



above Tl at the fixed point, i.e. 

+ 7^(2 + t^))^ = (k; 2t )^B t (<), (18) 
whose solution has the form: 



Hence, we can identify a^{j = 1) and ol^{t = 2) as equations (3) and (4). 

For the isotropic case, one has to consider only the momenta scale k 2 , corresponding 
to the quartic direction. The volume element in momenta space is A2. The (3 function, 
given by f3(u 2 ) = ~ e L 9l Q^ 2 , is different from the anisotropic one along the quartic direction. 
This implies that the isotropic critical behavior cannot be obtained from the anisotropic one 
whatsoever, for the coupling constants in both cases have different canonical dimensions. In 
the sequel, for sake of simplicity, one can replace (d— y) by | in the formulae for the critical 
exponents and 77^4 for the anisotropic behavior. Proceeding along the same lines, we 
obtain the scaling relations (0). 

From this analysis, it is easy to reproduce the critical exponents at one-loop level given 
in Ref. Jl| along the competing axes for both anisotropic and isotropic cases. Note that even 
though the fixed points are the same at one-loop level along the quartic and quadratic flow 
in the momenta scales for both cases, they do not have to be the same at higher loops. It 
follows naturally that a thorough description of the m-fold Lifshitz critical behavior actually 



needs four independent critical indices for the anisotropic case (in order to obtain all of them 
via scaling relations), reducing to two for the isotropic case. Their independence means that 
if we are willing to make suitable approximations to solve higher loop integrals, then they 
can be done independently in either subspace for the anisotropic case. 

We can check the validity of the hyperscaling relation (3) utilizing this reasoning. In the 
second paper of reference [0 the exponent vl2 and 7^2 were calculated within a new two-loop 
approximation. For d = 3,m = N = 1 (ANNNI model) the e^-expansion yielded 1/L2 = 0.73 
and T]i2 = 0.04. The exponent jL2 = 1-45 was obtained analytically from the e^-expansion 
together with the scaling relation. However, since = 1.5 is not a small parameter the 
neglected 0(e\) contribution could be relevant to the calculation of 7^2- In fact, had we 
replaced the numerical values of v L2 = 0.73 and r] L2 = 0.04 directly into the scaling law, we 
would have obtained 7^2 = 1-43. Therefore, whenever the expansion parameter > 1, one 
should use the numerical values of the exponents vl2 and 77^,2 obtained from the e^-expansion 
for fixed values of (m, N, d) in order to obtain the numerical values of the other exponents 
via scaling laws. The same reasoning applies when dealing with exponents associated to 
the competing axis, namely u L 4, f] L4: , etc. We can test this procedure by calculating other 
critical exponents perpendicular to the competing axis. By replacing the value of v^ 2 in 
(3), we find cti, 2 = 0.175, which is in remarkable agreement to the most recent Monte Carlo 



output ctL2 = 0.18 ± 0.02 Hl6[ . The susceptibility exponent calculated using this numerical 
method yielded 7^2 = 1-36 ± 0.03. On the other hand, from ([5c]), [3l 2 = 0.198, whereas that 
simulation resulted in /3 i2 = 0.238 ± 0.005, which is not as good as the result for the specific 
heat exponent, but is still reasonable. 



Indeed in flqj , the authors found the exponents along directions perpendicular to the 
competing axes and identified these exponents with a^, 7^ and After the work presented 
in this letter, we found appropriate to identify them with CK i2 > 7L2 and (3l2, respectively. 
Nevertheless, a curious fact takes place in the anisotropic case if in any given order in 
perturbation theory v L2 = 2v L4 and r/ L2 = \vla-, namely a L2 = a i4 = a Ll -y L2 = j L4: = -y L 
and /3l2 = @la = Pl- The scaling relations obtained in Ref. for the anisotropic case turn 

8 



out to reduce to those given here |TT]. This shows that the new results displayed here are 
consistent with precise numerical data. 

Previous RG treatments [1121 OI, introduced a dimensionful constant a in front of the 



first term in the Lagrangian presented in this work. In Ref. [12] another set of normalization 
conditions along with two symmetry points were chosen in order to renormalize the 1PI 
vertex parts. Those authors chose the symmetry point at nonvanishing quartic external 
momenta and zero quadratic external momenta for the vertex parts r^' 1 ' ) and T^' 2 \ 
In addition, they used two conditions on the derivative of the two-point vertex function at 
two independent external momenta scales. They argued that this choice of renormalization 
points would make the bare parameters and renormalization constants cr-independent. This 
reasoning implicitly takes into account the existence of only one fixed point, since it mixes two 

(2) 

different symmetry points at the same set of normalization conditions for T R , in conformity 
with the first approach based on a one-loop analysis [|]]. These two simultaneous conditions 
are responsible for the identifications 7^2 = 7l4 ; etc. The limitation of such renormalization 
group analysis is evident, since it fails to obtain the isotropic scaling relations. On the 
other hand, our analysis separates the two symmetry points in two independent sets of 
normalization conditions with independent renormalization group invariance in both sectors. 
Thanks to the condition 8q = fulfilled at the Lifshitz point, the first term of the above 
Lagrangian (quartic in the momenta) does not need to be multiplied by a dimensionful 
constant in order to make sense on dimensional grounds. Hence, a is not required at all in 
our new renormalization group approach, which permits freedom and simplicity to tackle 
the Lifshitz behavior in its full generality. 

The scaling relations presented for the isotropic critical behavior resembles the ones for 
the usual Ising-like system, the difference being the critical dimension in both cases. Except 
for the hyperscaling relation (which looks the same as the Ising-like one), all other scaling 
laws can be obtained from the usual (p 4 theory by making the replacement 2 — > 4. Besides, 
the isotropic behavior is completely independent of the anisotropic one. Then, the claim 



that one can treat the isotropic and anisotropic cases on the same RG grounds [0 is not 



consistent [|T^ . 

To summarize, we used a new renormalization group treatment for the m-fold Lifshitz 
critical behavior which constitutes a conceptual development in the comprehension of this 
sort of criticality. By using two different external momenta scales for the anisotropic behavior 
(perpendicular and parallel to the competition axes, respectively), we were able to derive 
two new sets of independent scaling laws, generalizing previous RG treatments. To our 
knowledge, we found for the first time new scaling relations for the isotropic critical behavior 
for m = d close to 8. A detailed description of the contents presented in this letter as well as a 
thorough description of all the critical exponents at least at two-loop order for the isotropic 
and anisotropic behavior will be given in a subsequent publication |IT |. The approach 



described here might be useful to treat Lifshitz points of generic character, when one allows 
arbitrary even momentum powers (greater than 2) in the Lagrangian (1), as well as other 
general field theories with higher-derivative appearing in other physical contexts ||. 
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